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Abstract 

Etching in semiconductor processing typically involves using halides because of the 
relatively fast rates. Bromine containing plasmas can generate high aspect ratio trenches, 
desirable for DRAM and MEMS applications, with relatively straight sidewalls. We present 
scanning electron microscope images for silicon-etched trenches in a HBr plasma. Using a 
feature profile simulation, we show that the removal yield parameter, or number of neutrals 
removed per incident ion due to all processes (sputtering, spontaneous desorption, etc.), dictates 
the profile shape. We find that the profile becomes pinched off when the removal yield is a 
constant, with a maximum aspect ratio (AR) of about 5 to 1 (depth to height). When the removal 
yield decreases with increasing ion angle, the etch rate increases at the comers and the trench 
bottom broadens. The profiles have ARs of over 9:1 for yields that vary with ion angle. To 
match the experimentally observed etched time of 250 s for an AR of 9: 1 with a trench width of 
0.135 pm, we find that the neutral flux must be 3.336 x 10 cm s . 


a) Electronic mail: hwang@dml.arc.nasa.gov 


I. Introduction 


Etching of silicon in the semiconductor industry is typically accomplished using halide 
containing plasmas, such as chlorine, fluorine, and bromine. Anisotropic profiles can be etched 
with ion-assisted etching using chlorine plasmas with fairly rapid etch rates. There are many 
advantages of using a relatively mature technology such as chlorine-based chemistry, such as fast 
etch rates, as well as having well-characterized material quantities and behavior (sputter yields, 1 
removal yields, 2 etc.). However, there are disadvantages to chlorine, such as poor selectivity to 
photoresist, as well as difficulty achieving high aspect ratio (AR) features. Other chemistries, 
such as bromine containing plasmas, have been demonstrated to etch highly anisotropic 
profiles 3 ’ 4 and are frequently used in commercial systems. Certain applications require extremely 
high AR devices, such as etching via holes, DRAM, and MEMS fabrication, which means that 
using solely chlorine-based chemistry is not viable. 

In bromine systems, the chemical mechanism of thermal etching has been studied, as 
well as recombination on the surface. 7 In spite of the fact that bromine has been used for more 
than fifteen years, its etch mechanism is not well understood and obtainin g material properties, 
such as etch yields as a function of incident ion angle and sticking coefficients of Br and Br 2 , is 
difficult. In general, bromine ion-assisted etch rates are seen to be lower than that for chlorine by 
a factor of about 1,6, 8 Although bromine etch yields are higher than for chlorine, it is thought 
that bromine adsorption onto the silicon surface has a considerably higher activation energy, 
which leads to a lower overall etch rate for bromine. 9 The total sputter yield for Br + has been 
measured to be higher than that for Cl + , although the chemical sputter yield for Cl is higher. 

In this article, we will describe the simulatioirwe used to calculate etched silicon trenches 
in bromine containing plasmas. We will show experimentally obtained high aspect ratio Si 
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trenches that were etched in HBr plasmas. We will show results from a feature profile evolution 
simulation that demonstrates the effect of ion removal yields on profile tapering and maximum 
obtainable AR. The yield magnitudes, as well as the functional dependence on ion impact angle, 
determine the width of the trench bottom. We will also make comparisons to the experimentally 
obtained SEMs. 


II. Experimental conditions 

The trenches are etched in Si with HBr plasma generated in a magnetically enhanced, 
capacitively coupled reactor with an applied frequency of 13.56 MHz. The inter-electrode gap is 
27 mm. The upper electrode has a quartz showerhead gas distribution plate to facilitate a 
uniform distribution of etching gases over the 200 mm diameter wafer. In addition, the reactor 
wall is outfitted with a quartz cylindrical shield for easy cleaning. The reactor itself is made of 
anodized aluminum. The reactor metal walls, as well as upper and lower electrodes, are at 60° C 
during processing. The wafer, however, reaches a different temperature depending on the 
process conditions, particularly RF power and reactor pressure. The cathode, or the lower 
electrode, also has a quartz focus ring 6 mm high, surrounding the wafer. The ring serves to 
reduce the electric field at the edge, resulting in improved etch uniformity. The wafer is held to 
the cathode by electrostatic force, the cathode being an electrostatic chuck with polyimide 
dielectric coating applied to the top surface of the cathode. An independent DC voltage is applied 
to the lower electrode to produce the desired electrostatic force. A magnetic field of 170 Gauss 
is applied by a set of rotating magnets, which facilitates a uniform field. The process conditions 
are 1300 W, 75 mTorr, and 175 seem HBr. An SEM sample of the etched trenches is shown in 
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Fig. 1, with a mask bevel of ±75°. The trenches are etched with high ARs of approximately 9:1, 
although the trenches taper at the bottom. The width at the bottom of the trench is roughly half 
of the width at the opening (surface), which is 0.135 pm wide. A close-up view of the trenches is 
provided in Fig. 2. The bottom of the trench is curved, which indicates that neutral shadowing 
decreases etching in the comers. This curvature suggests that the high AR of the trench limits 
the neutral flux at the bottom and lowers the neutral to ion flux ratio arriving from the plasma 
(r n /T + ~ 1). Also, since the trench is tapered, the ions are strongly directed toward the substrate 
with an anisotropic distribution function. The actual etch time was 4 minutes 10 sec (250 
seconds). 

m. Etch rate calculation 

Here we provide a brief summary of the approach we use to calculate trench etching. To 
describe the etch process, we assume a two step mechanism for etching of Si in HBr, similar to 
Cl etching of Si. This methodology is similar to that described earlier, 1 1 except here we 
distinguish between the two neutral species, Br and Br 2 . Passivation of silicon occurs due to Br 
and Br 2 reacting with Si to form SiBr x . Positive ions (HBr + , Br + ) impact the brominated surface 
and etching proceeds due to formation of volatile products. As in most commercial systems, 
obtaining plasma diagnostics was impossible due to reactor limitations; therefore, a simulation of 
the plasma conditions would lend few insights as there were too many unknowns to estimate and 
compare (species fluxes, electric fields, etc.). Instead, the species fluxes were calculated by 
estimating the distribution functions for the ions and neutrals. The distribution function for the 
neutrals was assumed to be Maxwellian. The ions were assumed to be Maxwellian in the bulk 
plasma and to fall collisionlessly through a sheath with voltage V s . Based on the tapered profiles 
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obtained in Figs. 1 and 2, the distribution functions for the ions have been selected to be strongly 
anisotropic and have the same analytic form given in Ref. 11, with eVJkT \ = 100. 

Using these estimated distribution functions, the fluxes must be calculated at each point 
in the 2-D (x, y ) computational domain around the trench. The directed velocity (of either ion or 
neutral) has the form |v|cos(0 -£), where |v| is the magnitude of the total velocity, 0 is the 

impacting angle from the normal, and £ is the angle of the surface from the normal. The 
impacting angle is defined as 


6 = tan 


-i 




\ Vy ) 


0 ) 


Then the directed flux from the plasma to a point on the surface is 


9 v 

r plasma = | Jjv|cos(0 -v)- /(jv|,0)-|v| d\v\ do . (2) 

0 

Note that Eq. 2 does not include re-emission fluxes from the trench surface. Both ion and neutral 

fluxes are calculated according to Eq. 2 for each species. The net ion flux is summed over all ion 

species and the neutral flux is calculated from the sum of the fluxes of Br and Br 2 . Re-emission 

fluxes for neutrals are neglected by setting the sticking coefficients equal to unity. This 

approximation represents the maximum possible etch rate by assuming that all neutrals will react 

with the surface. Ion re -emission fluxes are ignored, as microtrenching is generally not seen in 
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trenches etched in HBr plasmas compared to chlorine plasmas. 

In order to consider two separate neutral species, surface coverage parameters for both Br 
and Br 2 must be calculated. Experimental evidence has shown that bromine etch, like chlorine, 
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can be considered to occur in a monolayer, 4 ’ 13 hence using a surface site balance model is a valid 
assumption: 


<r Br = r Br s Br (1 -« Br -« Brz )- r + r Br a Br (3) 

at 

a,„ (l -a„ ~a Bli )- r ,Y B ,a Bri (4) 

Here, a is the surface coverage, a is the areal density, S is the sticking coefficient, T+ is the total 
ion flux, and Y is the removal yield, or number of atoms or molecules removed per incident ion. 
Note that the removal yield can be due to sputtering, spontaneous emission, etc. The surface 
coverages are found by assuming steady state and thus setting da/dt = 0. Each surface coverage 
parameter is then a function of removal yields and ion and neutral fluxes. Let 


d = r Br s Br r + r B + r Br s Br r + y Br +(r + ) 2 y Br y Bl 


(5) 


Then the surface coverages for Br and Br 2 are 
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Assuming that the etch yield of silicon (7 Sl ) is the same for all impacting ions, then the resulting 
etch rate is 
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( 8 ) 


ER=~ ^X(«B r +«Br 2 K,r + , 

Psi ions 

where psi is the mass density of silicon. The etch yield is dependent on impacting ion angle and 
has the functional form given in Ref. 1. 

Note that the properties that determine the uniqueness of a neutral species are its sticking 
coefficient, removal yield, and mass (which determine the distribution function). If S Br = S Br 2 
and T Br = T Br 2 , then Br and Br 2 are considered to be indistinguishable, etch wise. Furthermore, 
equations 6 and 7 then collapse to form the expressions for a and ER given in Ref. 1 1 for one 
effective neutral species (denoted n). 


a - 


r„ 


s„r n +y Bi t + 


(9) 


with an etch rate given by 


ER - 


■Y Si r + -a. 


( 10 ) 


IV. Simulation of Profile Evolution using Level Sets (SPELS) 

We have developed a feature profile evolution simulation, SPELS, which has been 
described earlier. 11 The level set equation is used to advance a higher order variable G(x, y, t), in 
which the plasma-surface interface occurs at G(x, y, to) = 0, or the zero level set. The plasma is 
defined where G(x, y, t) > 0 and the substrate by G(x, y, t) < 0. On the interface, the etch rate is 
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equal to the negative of the speed (5) of G. In the remaining (.x, y ) domain, G is updated by 
solving the level set equation, 

— + S|VG| = 0. (10) 

dt 1 1 

Therefore, the trench movement is indirectly tracked as G is advanced. The simulation continues 
for a user-specified time. 

It has been found that to obtain good resolution of the evolving feature, the mesh 
resolution typically requires about 80 x 80 cells per 1 x 1 normalized XIL, Y/L area. This would 
mean having 880 x 160 cells to model ARs of 9:1. Calculating the flux integrals over the entire 
mesh is computationally expensive and slows down the code significantly when the mesh is 
large. In order to run the simulation in a reasonable time, a narrow banding technique has been 
implemented. This technique is based upon the method derived by Adalsteinsson and Sethian 
and is outlined below. The level set equation is solved only in a “band” around the zero level set, 
thus reducing the number of “active” cells. G is set to a minimum (negative) value beyond the 
edge of the band in the solid, and a maximum value in the plasma above the band edge, as shown 
in Fig. 3. The band is reconstructed whenever the zero level set gets close to the edge of the 
band. For the following cases, the band is approximately 22 cells wide (1 1 cells above and 
below the zero level set), and the band is reconstructed whenever the zero level set gets within 6 
cells of the edge of the band. Using this method, the computational speed-up is a factor of 5 or 
more. 

Obtaining most physical quantities for Br etching of silicon is difficult, as many of the 
parameters are unknown for HBr interactions with Si. Asa first attempt, we have assumed 
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material parameters that are similar to chlorine, given in Ref. 1 1, unless specified otherwise. We 
also assume that the form of the silicon etch yield versus ion impact angle for Br and Br 2 is 
similar to that of chlorine given in Ref. 1, although calculations involving the relative etch yield 
for HBr is considered later in this article. In the following cases, we do not consider diffusive re- 
emission of neutrals and instead set the sticking coefficients of Br and Br 2 to unity. The 
normalized length scale is L = 0. 135 pm. The amount of recombination of bromine on the 
reactor surfaces is unknown and we assume that the ratio of neutral fluxes from the plasma Tg^/ 
Ter = 3 and the ratio of the ion fluxes is 1^+/ THBr+ = 3. The ratio of the total neutral flux to the 
total ion flux is unity (T n / T+ = 1). The hard mask is assumed to be a height of 0.2 L, which 
corresponds to about ±78° angle allowance at the center of the trench, which is similar to the 
experimental mask bevel allowance. 


V. Results and Discussion 

The following results were obtained using the formulation given in Eqs. 6-8 for the 
surface coverages and etch rates. In Fig. 4(a), Ter = Ybi2 — 1 (meaning Br and Br 2 have identical 
properties), leading to a maximum aspect ratio of 5:1 for the trench. The neutrals become 
shadowed inside the trench and are less able to brominate the surface, particularly in the comers. 
The etch rate subsequently slows down at the comers and the trench pinches off. Increasing the 
removal yields to Ter = Ybt 2 = 3, as shown in Fig. 4(b), broadens the trench slightly but the 
maximum AR is still not substantially greater than 5:1. 

The trenches shown in Fig. 4 assume that the removal yields Fe r and Yq t 2 are constant 
with impacting ion angle. This assumption may not be correct, as the removal yield encompasses 
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all methods of removing the brominated surface species, including physical sputtering. For 
example, ion assisted etch yields have been shown to be a strong function of ion angle (Ref. 1), 
and are included in calculating F' r+. The trenches shown in Fig. 5 employ F Br and F B r 2 which 
vary with impacting ion angle in the same way that Fsi varies. In Fig. 5(a), the maximum 
removal yield F Br = F Br 2 = 1, and in Fig. 5(b), the maximum yield is F Br = F Br2 = 3, corresponding 
to cases in Figs. 4(a) and (b), respectively. In these cases, using a yield that varies with ion angle 
leads to increased surface coverages in the comers. The higher surface coverages create broader 
trenches and thus generate higher aspect ratios. The trench in Fig. 5(a) is more square than the 
pinched off trench in Fig. 4(a), although there is some bowing of the sidewalls. Increasing the 
removal yield in Fig. 5(b) leads to a flared trench, which has a severe amount of undercutting. 
Both trenches have rounded bottoms, which indicates that the surface coverages vary along the 
bottom due to neutral shadowing. Higher aspect ratios will increase the amount of shadowing. 

Allowing an angular dependent removal yield essentially scales the total ion flux and 
modifies the surface coverages in the following way. Both a Br and a Br2 can be calculated using 
Eqs. 6 and 7, but instead have the following expression instead of Eq. 5: 


d — r Br 5 Br r + F B 


T ^Br, ^Br,r+^Br 


+ 0(rJ 


Y Y 

1 Br^Br, 


( 11 ) 


where 0 < /? < 1 . Hence surface coverages that depend on ion angle will always be equal to or 
larger than the coverages using constant removal yields. The amount of “flare” in the trench can 
be dictated by varying the dependency of F Br and F Br2 on ion angle. For example, one can select 
an artificially strong dependence on ion angle such as 


Y„ =F o exp(-|0-£|). 


( 12 ) 
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This extreme dependence leads to undercutting as shown in Fig. 6. Allowing the removal yields 
of both Br and Br 2 to vary according to Eq. 12 increases the surface coverages in the comers and 
causes severe undercutting. In contrast, the comers become shadowed and do not etch when the 
removal yield is set to a constant value (F Br = F Br 2 = 1). This constant value yield case is the 
same as shown in Fig. 4(a); these two cases represent the two extremes for the functional 
dependence of the removal yield on impact angle. The example shown in Fig. 6 and given in Eq. 
12 is for demonstration purposes only, as the removal yield is expected to depend on ion angle in 
a similar fashion to the ion-enhanced sputter yield. When the removal yield is constant, the 
neutrals do not sufficiently brominate the comers and thus the trench pinches off, creating a low 
AR trench. 

The following cases use a removal yield for Br and Br 2 that varies with ion impact angle. 
The yields given are for the maximum values, and the functional form varies in the same manner 
as does the etch yield Fsi. Including an angular dependence of the removal yields can cause a 
flared shape, as seen in Figs. 5(b) and 6. However, the amount of tapering in the trench can also 
be adjusted by scaling the maximum value of F Br . For large removal yields (F Br = F Br 2 = 3) as 
seen in Fig. 7 (a), the trench is undercut much in the same way that as is the profile in Fig. 6. In 
other words, assuming the removal of Br and Br 2 on the surface is less preferential for ions that 
impact at grazing angles causes the trench to etch more in the comers and thus leads to the flared 
shape. The tapered trench shown in Fig. 7(b) results from a low removal yield, or F Br = y Br 2 = 
0.4, which essentially sets the ion-assisted etch yield equal to the removal yield (no spontaneous 
desorption of neutrals). 
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The above cases all assumed that both Br and Br 2 had identical sticking coefficients and 
removal yields. The trenches shown in Figs. 8(a) and (b) assume different removal yields for Br 
and Br 2 . For both cases, r Br 2 /r Br = 3, meaning a low degree of dissociation. The trench in Fig. 
8(a) has Y B , = 0.4 and F Br2 = 3, which leads to a flared trench as the majority of the neutral flux is 
Br 2 which has a high removal yield. The tapered trench in Fig. 8(b) has the opposite yields, or 
Y Bi = 3 and F Br2 = 0.4. Again, the dominant neutral flux and its removal yield will dictate the 
shape of the trench. 

Vitale et al . 15 measured relative etch yields for HBr and Br 2 and found that the yields 
have a weaker dependence on ion impact angle than does the yield for Cl 2 , which was originally 
measured in Ref. 1. These differences are shown in Fig. 9. The etch yield for Cl is constant for 
small angles of incidence, but then has a sharp drop-off for large angles. In contrast, the etch 
yield for HBr has a slow decrease with increasing ion angle. However, it should be noted that the 
authors indicate that the experimental error is significant for larger angles (approximately ±53% 
for HBr compared to ±17% for Cl 2 ). The angular dependencies for etch yields for chlorine and 
HBr fall within each other’s error ranges, making it difficult to determine if the overall shapes are 
correct. In any case, for high aspect ratio cases the overall angular dependency of the etch yield 
may not be significant in determining profile shapes. The case shown in Fig. 10(a) uses Vitale’s 
etch yield and has the same overall profile shape as does the case in Fig. 10(b), which uses 
Chang’s yield for chlorine. The magnitude of the etch yield in Fig. 10(a) is twice that in Fig. 
10(b); however, the tapered trench shape is indicative of etching that is limited by neutral 
shadowing. Both cases assume that the removal yield has the same angular dependence as their 
respective etch yields. Because ion reflections are not considered, the profile shape is not 


11 


strongly dependent on the angular dependence of the etch yield, although the overall etch rate 
does increase. 

Although the overall AR obtained from the model can be matched to the experimental 
SEMs, the profile shape differs for the bottom portion of the trench. This difference can be seen 
clearly in Fig. 11, where the model results from Fig. 7(b) are superimposed on the SEM. The 
etched trench has a gradual tapering, whereas the predicted profile is broader and suddenly tapers 
at the bottom. Using an ion-assisted yield of Fsi = 0.4, the total neutral flux required to match the 
experimental etch time is 3.336 x 10 17 cm 2 s *. Several possibilities exist as to why the profile 
shapes do not match exactly. Other phenomena, which were not considered in our simulation, 
may contribute to the overall profile shape. For example, it is possible that ion reflection may 
contribute to etching. Vyvoda et al. 16 theorized that in HBr plasmas, ions reflect off of the 
silicon in a diffuse distribution. Although their model overpredicted the amount of etching in the 
comers for a trench with AR -3.3:1, they did find that their results were sensitive to the assumed 
reflected ion distribution. Thus, reflected ions contributing to the etch process, particularly in the 
comers, may account for the difference in the actual profile versus the calculated one. Other 
secondary processes, such as H-atom enhancement of the etch rate due to coadsorbtion, lattice 
penetration, and Si backbond breakage, 15 may also contribute to the broader trench bottom 
profiles. Another possibility is that our calculation is 2-D in nature, and solving for a full 3-D 
calculation of the etch may improve the profile shape. 

VI. Concluding Remarks 

Trenches with aspect ratios of 9: 1 were etched in HBr plasmas. Incorporating a narrow 
band technique allowed for calculations of large geometries that would otherwise not be possible. 
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This is particularly important for high aspect ratio etching, where trenches can exceed ARs of 
20:1. For the calculated profiles, the removal yield of a neutral on silicon was found to be 
important for it can dictate the profile shape. Higher yields lead to broader trench bottoms (and 
ultimately to higher aspect ratios). If the yield has a functional dependence on ion impact angle 
in the same manner as the ion assisted etch yield, the trench bottom can become significantly 
wide and become flared. Although the overall AR can be matched to experimental values, many 
parameters (sticking coefficients, species fluxes, etc.) are unknown and were estimated, thus the 
predicted profile shapes are different. This also emphasizes the need to integrate the feature 
evolution model to a reactor model self-consistently in order to obtain various fluxes; however, 
the current limitation is the lack of knowledge on collisional cross sections for electron and ion 
impact reactions and details of surface chemistry. 
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Figure Captions 


Fig. 1. SEM of etched trenches of silicon obtained in HBr plasma reactor. Aspect ratio is 9: 1 per 
trench with etch opening of 0.135 pm. Total etched time is 250 seconds. 

Fig. 2. Close-up view of trenches. Trenches have gradually tapered sidewalls and curved 
bottoms, which indicate neutral shadowing. Neutrals are shadowed both due to the high AR and 
because of low neutral flux relative to the ion flux from the plasma. 

Fig. 3. Sample of contours from SPELS which define the location of the plasma-substrate 
interface. The narrow band technique only defines G “near” the front, (a) For plasma locations 
far away from the surface, G = G max . The substrate far from the surface is defined by G = G m i n - 
(b) A close-up view of the contours shows that within the band, G > 0 defines the plasma region, 
G = 0 defines the surface, and G < 0 indicates the substrate. 

Fig. 4. Profiles generated from SPELS. Pinched off trenches with maximum AR 5:1, with 
constant values of removal yields. Total ratio of neutral flux to ion flux is unity. Neutrals are 
shadowed in the comers which prevent obtaining higher ARs. (a) TBr = Tbf 2 = 1, and (b) Ter = 
F B r2 = 3. 

Fig. 5. Same conditions as in Fig. 4, except removal yields are functions of ion impact angle. 
Less removal of Br and Br 2 in the comers of the surface allow for deeper etching, (a) FBr = TBr 2 
= 1, and (b) F Br = F Br 2 = 3. 

Fig. 6. Difference between constant removal yield function and one that varies with ion impact 
angle. The dotted line is for the same case in Fig. 4(a), and the solid line has the extreme 
functional dependence of removal yield on ion impact angle given in Eq. 12. 

Fig. 7. Tapered and flared trenches with non-constant removal yields. Both removal yields have 
functional forms similar to etch yield of Si with Cl, as given in Ref. 1. (a) Tar = FB r 2 = 3 and (b) 
F Br =F B r2 = 0.4. 


Fig. 8. Profiles assuming different removal yields for Br and Br 2 . For both cases, Tb^/T B r = 3. 

(a) y Br = 0.4 and F B r2 = 3. (b) Y Bi = 3 and F Br 2 = 0.4. The dominant neutral flux and its removal 
yield define the shape of the trench. 

Fig. 9. Experimentally observed silicon etch yields as a function of impacting ion angle, taken 
from Refs. 1 and 15. The etch yield in HBr plasmas is weakly dependent on ion angle, whereas 
the etch yield in chlorine plasmas demonstrates a sharp decrease with large angles. Error 
estimates for the HBr etch yield at large angles is ±53%, and ±17% for chlorine for mid-range 
angles. 

Fig. 10. Comparison of profiles using etch yields shown in Fig. 10. (a) Etch yield using angular 
dependency for HBr with F Br = T Br 2 = 0.2. (b) Etch yield using chlorine angular dependency, 
with F B r = Tbf 2 = 0.4 (same as in Fig. 7(b)). 

Fig. 11. Superimposed model results from Fig. 7(b) on the SEM from Fig. 1. The aspect ratio 
can be matched but the profile shapers differ, largely in part due to estimated yields as well as ion 
and neutral fluxes. 
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